When can Fokker-Planck Equation describe anomalous or chaotic transport ? 
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The Fokker-Planck Equation, applied to transport processes in fusion plasmas, can model several 
anomalous features, including uphill transport, scaling of confinement time with system size, and 
convective propagation of externally induced perturbations. It can be justified for generic particle 
transport provided that there is enough randomness in the Hamiltonian describing the dynamics. 
Then, except for 1 degree-of-freedom, the two transport coefficients are largely independent. De- 
pending on the statistics of interest, the same dynamical system may be found diffusive or dominated 
by its Levy flights. 
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The Fokker-Planck Equation (FPE) is a basic model 
for the description of transport processes in several sci- 
entific fields. In one dimension it reads 



dMx.t) = -d., {V{x)n) + dl {D{x)n) 



(1) 



where n is the density for a generic scalar quantity, V the 
dynamic friction, and D the diffusion coefficient. In par- 
ticular, FPE backs up the diffusion-convection picture of 
anomalous transport in magnetized thermonuclear fusion 
plasmas, where V is often referred to as a pinch veloc- 
ity. This transport stems from turbulence triggered by 
density and temperature gradients through convective in- 
stabilities. Though very popular, the drift-diffusive pic- 
ture underlying FPE breaks down in some cases. This 
was proved, e.g., for the transport of tracer particles 
suddenly released in pressure-gradient-driven turbulence, 
which exhibits strongly non gaussian features fl|. This 
fact triggered a series of studies where transport was de- 
scribed in terms of continuous random walks with Levy 
jumps, and of fractional diffusion models (see and ref- 
erences therein). This sets the issue: when is FPE rele- 
vant for anomalous or chaotic transport, when is it not? 

Our discussion has many facets, but goes through the 
following main steps : (i) FPE with a source term can 
model phenomena commonly labeled as "anomalous" in 
fusion plasmas: uphill transport Q, anomalous scaling 
of confinement time with system size and non dif- 
fusive propagation of externally induced perturbations 
(ii) If the source is narrower than the mean ran- 
dom step of the true dynamics, FPE fails, while the 
correct Chapman-Kolmogorov equation reveals the ex- 
istence of spatial features of the spreading quantity that 
are not related to the transport coefficients. {Hi) For 
general 1 degree-of-freedom (dof) Hamiltonian dynamics, 
the constraint V — dD/dx holds, as originally derived by 
Landau. This constraint vanishes for higher dimensional 
dynamics, incorporating for instance particle-turbulence 
self-consistency, or the full 3-dimcnsional motion of par- 
ticles, (iv) FPE can be justified for particle transport 
provided that there is enough randomness in the Hamil- 



tonian describing the dynamics: e.g., when it includes 
many waves with random phases. This may work even 
whenever the dynamics of individual particles exhibit 
strong trapping motion. However chaos is not enough to 
justify FPE. (v) The diffusion may be quasilinear (QL) or 
not, depending on the so-called Kubo number K, which 
scales like the correlation time @, 01 • For K small, QL 
diffusion is due to the locality in velocity of wave-particle 
resonance. Large K diffusion is due to the locality of 
trapping in phase-space. 

Our discussion aims at providing insight in transport 
mechanisms, rationales for confinement scaling laws, and 
tools for experimental data analysis. It involves several 
building blocks: some are provided by the available liter- 
ature, while a few are elaborated here. Though written 
for the fusion community, a large part of our discussion is 
of direct relevance to chaotic transport, and in particular 
to Lagrangian dynamics in incompressible turbulence. 

What FPE can do. The case where V and D are con- 
stant in Eq. ([T]) is often considered in fusion data analy- 
sis. An initially Dirac-like perturbation travels with ve- 
locity y, while diffusing with coefficient D. For large 
enough systems, L ^ D/V, the perturbation actually 
travels ballistically. Thus FPE can model the non dif- 
fusive propagation of induced perturbations in a fusion 
machine ,5]. The scaling of confinement time Tconf with 
system size L, goes from diffusive (rconf « i^, L ^ D/V) 
to ballistic (rconf oc i\ L > D/V). Thus FPE may ac- 
count for anomalous scaling of confinement time with 
system size. It is even more so when V and D depend 
on cc or when their possible dependence on spatial 
gradients driving turbulence is accounted for. Whatever 
V and D be, FPE may be written as dtn — —dxT, with 
the fiux r = Vn — dx{Dn) . For D constant, if there 
are no sources in the central part of the plasma, for sym- 
metry reasons F = 0, which yields dxn/n — V/D: if 
sign{x)V < 0, where x — corresponds to the plasma 
center, the probability piles up toward small |a;|'s, leading 
to an equilibrium distribution of "uphill transport" type 
3]. However, the piling up may result as well from V = 0, 
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and D{x) growing with \x\ since n{x) oc 1/D{x). This is 
a caveat for data analysis: a broad family of {V, D) pro- 
files can model the same experimental data. Let us notice 
that the ability of FPE to describe anomalous transport 
is shared with models giving non gaussian features (see 
and references therein). 

Transport in presence of a source. The classical 
derivation of FPE from the Chapman-Kolmogorov equa- 
tion 

dtn{x,t)= f P(x,x')^^^dx'-^^ + S{x), (2) 
J t{x') t{x) 

where S{x) is a localized source term and t{x) is a wait- 
ing time, assumed to be 1 for the moment, assumes the 
typical width ctp of P in a; — x' to be the smallest spa- 
tial scale of interest. Therefore the width as of S must 
be much larger than ap for FPE to describe correctly 
n{x, t) within the source domain; otherwise the source 
pumps strong gradients in n{x, t), and the Taylor expan- 
sion leading to FPE breaks down close to it. Assume 
P{x,x') to be a function of x — x' only. Then Eq. ^ 
shows that for scales smaller than ap, the Fourier trans- 
form of a stationary n{x) is almost equal to that of S{x). 
These are the scales relevant to describe the profile close 
to the source. Therefore, n{x) displays a bump similar 
to that of S{x) that is not due to specific properties of 
the FPE transport coefficients: a further caveat for data 
analysis. 

Paradigm model. Consider tracer transport in 2- 
dimensional (2D) incompressible turbulence, or transport 
in magnetized plasmas induced by electrostatic turbu- 
lence in the guiding center approximation. A paradigm 
for such a transport is the equation of motion 

^--Vcf>(X,i) xx3^v (3) 

with X = {xi,X2) the particle position (perpendicular to 
the magnetic field, for the plasma case) and $ the flow 
function or the appropriately normalized electrostatic po- 
tential. Model ^ applies to transport due to magnetic 
chaos as well. Indeed, for an almost straight magnetic 
field, $ may be replaced by <&' = $ — where both 

the electrostatic potential $ and the parallel vector po- 
tential are computed at the guiding center position, 
and f|| is the parallel velocity [6|]. Therefore Eq. ^ 
describes the canonical equations for the guiding center 
dynamics in a mixture of electrostatic and magnetostatic 
fields, ruled by Hamiltonian $(X, t), with conjugate vari- 
ables xi and X2. 

Constraint on FPE due to the dimensionality. Con- 
sider model ([3]) where $ has a bounded support, or is 
periodic in the Xi's with an elementary periodicity cell 
R: motion is, or can be considered as, located within a 
bounded domain R of phase space. Because of conserva- 
tion of the area of a phase space element during motion, 
an initially uniform particle density n in R must remain 



uniform for later times. Now, assume FPE describes the 
evolution of n{x) {x one of the two conjugate variables). 
Since n is stationary, the flux F must be a constant in 
i?, and must vanish, since it does on the boundary of 
R. Since F — Vn — dx{Dn), this requires V = dD/dx. 
This is the fundamental reason for this constraint, first 
derived by Landau Q for a stochastic but non chaotic 
Hamiltonian dynamics, since his derivation uses Taylor 
expansions of the orbit in time. The Landau constraint 
(LC) was recovered more recently for the chaotic motion 
of particles in a prescribed set of Langmuir waves Q . If 
V — dD/dx, then F = —Ddn/dx. In absence of sources 
in the plasma core, F = implies dn/dx in this do- 
main, which rules out uphill transport. 
Adding more dimensions to the previous Hamiltonian dy- 
namics (for instance the parallel motion in addition to 
the ExB drift) leads in general to a breakdown of LC. 
Indeed a Kolmogorov-Arnold-Moser torus is no longer 
able to separate the chaotic domain into disconnected 
sets (Arnold diffusion) , and therefore there are no longer 
boundaries F must vanish onto. For instance, LC no 
longer holds for the self-consistent motion of particles 
in a set of Langmuir waves (here x is the particle ve- 
locity): V is equal to dD/dx plus a drag force due to 
the spontaneous emission of waves by particles p. The 
analogy of Langmuir wave-electron interaction with the 
toroidal Alfven eigenmode-fast ion one, where x is the 
radial position *lCr|, shows that a similar effect may hold 
for a fusion machine. This is a particular instance re- 
vealing that the true Hamiltonian dynamics of particles 
in fusion machines has more than 1 dof. In general LC 
cannot hold then: uphill transport and a central finite 
density gradient are possible. 

Derivation of FPE from Hamiltonian dynamics for 
particles. We consider model ^ where $ is a statisti- 
cally stationary, spatially homogeneous, isotropic zero- 
mean-value stochastic potential with typical amplitude 
<f>Oj a-iid a given two-point, two-time correlation func- 
tion, with correlation time Tc and correlation length Ac, 
as considered in Refs. [1, [llj. Some results from these 
references are important for justifying FPE in this case. 
Potential $ drives the particles motion by setting the 
instantaneous value of their velocity, whose typical am- 
plitude is $o/Ac. The system chooses among two types of 
diffusive chaotic dynamics according to the value of the 
Kubo number K = ^oTc/A^. The rationale for this is 
the following. If the potential is static (tc infinite), par- 
ticles are trapped into potential wells and hills, and pos- 
sibly make long flights along "roads" crossing the whole 
chaotic domain. These various domains are separated by 
separatrices joining nearby hyperbolic points. If Tc is fi- 
nite, but large {K ^ 1), the potential topography slowly 
changes, and the dynamics evolves quasi-adiabatically. 
Since phase space area inside the instantaneous closed or- 
bits must be adiabatically preserved, and since the area 
of the various domains defined by the separatrix array 
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FIG. 1: (a) 3> 1: as the net of intersecting separatri- 
ces evolves in time (from solid to dashed curves), the lower 
domain enlarges at the expense of the upper one. Hence a 
trajectory may go from the upper to the lower one (thick to 
dashed curve), (b) K <^ 1: jumps among small trapping arcs 
in a quickly evolving potential topography. Solid and dashed 
contours stand for potential hills and wells. 



fluctuates a lot, orbits must cross the instantaneous sep- 
aratrices, and jump this way from one domain to the next 
one (Fig. [T{a)) . In the absence of "roads", which is al- 
most the case for a gaussian spatial correlation function 
of the potential [llj, these crossings produce a random 
walk with step Ac and waiting time Tc; the correspond- 
ing diffusion coefficient is I? « A^/tc. The presence of 
"roads" modifies this estimate and brings some depen- 
dence upon K ^1^. As a result, for ivT ^ 1 diffusion is 
justified by locality of trapping in phase-space. 

When K I, the particles typically run only along a 
small arc of length ^qTc/ Ac of the trapped orbits of the in- 
stantaneous potential during a correlation time. During 
the next correlation time they perform a similar motion 
in a potential completely uncorrelated with the previous 
one (Fig. [Ub)). These uncorrelated random steps yield 
a 2D Brownian motion with a QL diffusion coefficient 
D w ^Itc/XI- These two limit cases in K show diffusion 
is a quite general behavior of particle transport, even 
whenever structures are visible in the electrostatic po- 
tential, or whenever non gaussian behaviour is obtained 
for a more limited statistics, as in Ref. [ij whose dynam- 
ics, except for isotropy, may be thought as one realization 
of that in Ref. 
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The simple preceding reasoning strongly bears on the 
stochastic properties of the potential, and not on the 
chaotic features of particle dynamics. A more rigorous 
picture for the diffusion process of the K <^ 1 regime 
can be given, which incorporates chaos as an essential 
ingredient, but exhibits the paramount importance of 
potential randomness. This picture is a translation for 
dynamics ([3]) of that described in Refs. for the 

dynamics of an electron in a set of Langmuir waves. To 
this end, we consider $ as a sum of propagating modes 



eaicw cos (k • X — oji -I- i^k.Lij), where the (/3k,uj's are uni- 
formly distributed random phases, and the spectrum is 
isotropic in k; e is put for scaling purposes. If the parti- 
cle does not resonate with the modes in ([3]), its dynamics 
may be described by perturbation theory in s. In the op- 
posite case, at any given time to where the particle has 
velocity v(to), some modes are resonant with the parti- 
cle, and their action cannot be described in a perturba- 
tive way, but some are not, and still act perturbatively. 
These two classes can be defined according to their reso- 
nance mismatch with the particle p(to) = |k ■ v(to) — 
0, 12, 13 1 . Consequently, the instantaneous particle mo- 
tion splits into a perturbative, and thus non chaotic part, 
and into a non perturbative chaotic part due to the set 
S'(v(to)) of modes which are resonant enough with the 
particle at time Iq: p{tQ) less than some threshold po e. 
This set evolves with time, according to the instanta- 
neous value of v{t). A suitable choice of po/e enables to 
incorporate the set of modes driving the chaotic dynamics 
of the particle over a finite interval [tQ—At, to+At], where 
At is of the order of the time Tsproad ~ X^/D. Over this 
time interval the chaotic particle dynamics is ruled by a 
reduced Hamiltonian ^(^(X, i), which is 3>(X, t) with the 
summation restricted over the modes inside set S'(v(to)). 
If the dynamics is chaotic, v(t) changes a lot its direc- 
tion during its motion, which brings disconnected sets 
S'(v(i„)), for different times tn,n — 0,1,.... The ran- 
domness of the (^k.w's implies that the dynamics ruled 
by each of the <I>t^(X,t)'s are statistically independent, 
which justifies a central limit argument for their cooper- 
ative contribution to the particle motion: the dynamics 
is indeed diffusive. 

The above argument uses the locality of wave-particle 
resonance in phase velocity. It also holds to prove the 
diffusive behavior for peaked frequency spectrum, or for 
a single typical outcome of the random phases, if the set 
of the initial particle velocities vq is spread enough for 
them to be acted upon by a large number of disconnected 
S'(vo)'s. This explains the diffusive behavior found in 14 1 
by averaging over initial particle positions. It should be 
stressed that the above locality rationale depends essen- 
tially on the random phases, and that chaos, through the 
breakup of KAM tori, just brings in the ability of the 
motion to be ruled successively by uncorrelated dynam- 
ics. In the absence of random phases no diffusivity can 
be derived by the above reasoning. Indeed the diffusion 
picture was clearly shown to break down inUj, and for 
electron dynamics due to Langmuir waves [l3|. 
The existence of random phases can also be used to de- 
rive a rigorous QL estimate for D. This derivation is 
analogous to that for the chaotic dynamics due to Lang- 
muir waves for A' ^ 1 9]. Its central argument is that 
the dynamics depends slightly on any two phases during 
a time much larger than Tspicad, which defines the time 
of strong sensiti vity of the dynamics on the whole set of 
random phases [I4I . These successive two random-phase 
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arguments for the K <^ 1 case, assumed for simplicity 
that all phases were uncorrelated, but some correlation 
may be accommodated. We stress they do not use at all 
any loss of memory due to chaotic motion: indeed differ- 
entiable chaotic Hamiltonian dynamics is not hyperbolic. 
Beyond simple Hamiltonian models. More recently 



Vlad et al. [15| introduced a spatial inhomogeneity in 



model ([3]), such that its r.h.s. is multiphed by a grow- 
ing function of xi. This was meant as a modeling of the 
increase of the magnetic field toward the main axis of 
a fusion machine, but makes the dynamics non Hamil- 
tonian stricto sensu. On top of the previous diffusive 
behavior (which becomes anisotropic), the new inhomo- 
geneity brings a "radial" drift velocity V along xi due 
to the chaotic motion, which corresponds to the dynamic 
friction of FPE. The sign of V depends on K. If if < 1, 
since the velocity increases toward larger xi's, the dis- 
placement during a correlation time is larger toward the 
exterior than toward the interior, which brings an outgo- 
ing drift. If » 1, the trapped particles are slower in 
the inner part of their orbit, which increases their prob- 
ability to be there with respect to that to be in the outer 
part: this brings an ingoing drift. Since D now grows 
with xi,V — dD/dxi is impossible for K ^ I. 
The previous discussion can be extended to higher di- 
mensional systems quite naturally, as do the above ar- 
guments of locality, or the argument of the weak effect 
of two phases on the dynamics. The diffusive aspect of 
higher dof chaotic systems is largely documented as well 
(see [l^ and references therein). In conclusion we may 
state that depending on the statistics of interest, the same 
dynamical system may be found diffusive or dominated by 
its Levy flights. In particular averaging over many ran- 
dom phases may lead to diffusion, while averaging only 
over a limited set of initial conditions for the particles 
may lead to the Levy flight picture. Since fusion exper- 
iments generally average over many plasma realizations, 
and global confinement scaling laws even more so, FPE 
is a highly relevant tool for this field. It should be noted 
that in a magnetized toroidal plasma, density n describ- 
ing particle transport is the true particle density divided 
by a growing function t{x) of the local magnetic field 
(see [l7| and references therein). This brings a slanting 
of the density profile toward the outer part of the torus 
that has nothing to do with a turbulent transport phe- 
nomenon, in contrast with that in Ref. [isj . and which 
corresponds exactly to a waiting time t{x) in Eq. ([2]). 
This is one more caveat for data analysis. 
We found that for ^ 1, FPE with a QL diffusion coef- 
ficient is justified by chaotic Hamiltonian dynamics with 
random phases, even though structures exist in phase 
space for one realization of the phases. This, and the 
fact that the QL diffusive modeling of transport is quite 
efficient [3], suggest that K is small in magnetic fusion 
turbulence. There are reasons for K not to be large. 
First, the usual estimate in fluid mechanics of the cor- 



relation time as the eddy turn-over time yields K = 1 
19]. Furthermore strong turbulence theory predicts that 
K is at most of order 1 (see Eq. 4.34 of (l^l)- How- 
ever the issue of the typical value of K is still unsettled, 
since an analysis of fluctuation data in the TEXT toka- 
mak indicated K of order 1, and a. poor agreement of QL 
estimates for impurity transport 21|. It would be inter- 
esting to systematically compute K from experimental 
or numerical data. In rotating or stratified fluid turbu- 
lence a weak effect of structures holds as well: cigar-like 
or pancake-like structures are present, but turbulent dif- 
fusion is correctly modeled by assuming random phases 
for the Fourier components of the turbulent fluid [24I . 
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